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, Abstract 



Q"^ I In a schematic model with equidistant fourfold degenerate single- nucleon levels, 

a conventional isovector pairing force and a symmetry force, the RPA correlation 
energy rises almost linearly with the isospin T, thus producing a Wigner term in 
Q"^ ' accordance with the empirical proportionality of the symmetry energy to T{T + 1). 
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Nearly symmetric nuclei have an extra binding, the so-called Wigner energy, 
that is not described by the quadratic symmetry term in the semi-empirical 
^ I mass formula [1]. This is explained in various ways in the literature. Counting 

' the even bonds among supermultiplet degenerate nucleon orbitals, Wigner 

estimates that the isospin-dependent part of the interaction energy in the 
ground state of a doubly even nucleus is proportional to T(T -|- 4), where T 
is the isospin [2]. Talmi proves that for seniority-conserving forces such as the 
pairing force acting in a single j-shell, this part of the interaction energy is 
proportional to T(T + 1) [3], and Bohr and Mottelson point out that this 
isospin-dependence arises in general from an interaction proportional to the 
scalar product of the isospins of the interacting nucleons [4]. A symmetry en- 
ergy proportional to T(T + 1) also fits the empirical masses well [5]. Myers 
and Swiatecki attribute the extra binding of nearly symmetric nuclei to the 
interaction of neutrons and protons in overlapping orbitals [6]. Shell model 
calculations with realistic forces are succesful in reproducing the measured 
binding energies [7], whereas with Skyrme forces, no Wigner term appears in 
Hartre-Fock-Bogolyubov calculations and only a small one in Hartree-Fock 
calculations [8]. However, by invoking a particular isoscalar pairing force that 
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breaks geometric symmetries, Satula and Wyss obtain a significant Wigner en- 
ergy in approximately number-projected Bogolyubov calculations [9]. In finite- 
temperature Bogulyubov calculations with a Yamaguchi force, Ropke et al. 
get at low temperature in the local density approximation a contribution from 
neutron-proton pairing to the binding energies in the A = 40 isobaric chain 
with a maximum aX, N — Z — —1 and a quite irregular dependence on the 
asymmetry [10]. 

The RPA is the leading order correction to the Hartree-Fock-Bogolyubov ap- 
proximation. I have therefore calculated the RPA correlation energy from the 
schematic Hamiltonian 

H^Ho-GP^-P+ //o = E efc«L%.r ■ 

kar 

In this expression, the index kar labels orthonormal nuclcon orbitals, and 0^0-7- 
are the corresponding annihilation operators, ka takes the values k and k so 
that the orbital k is obtained from the orbital k by time-reversal, and r is 
'n' for a neutron orbital or 'p' for a proton orbital. P is the pair annihilation 
isovector, and T denotes the total isospin. The former has the coordinates 

k k 

The single- nucleon energy takes equidistant values separated by 77, and 
G and k are coupling constants. 

To describe states with a given number = J2kaT ^kar^kar of valence nucleons 
and a given isospin, I employ the Routhian 

R = H -XA^- i^T,. 

It is just for convenience that is chosen here as the isospin-coordinate to be 
constrained. Since H is isobarically invariant, one could equivalently constrain 
the projection of T on any axis in isospace. Following Marshalek [11] I base 
the RPA on the Hartree-Bogolyubov (not Fock) self-consistent state derived 
from R. This is the Bogolyubov vacuum that minimizes £'0 — X{Ay) — iJ,{Tz), 
where 

£^0 = W-^ + f(T)^ A = -G(P). 

At the minimum one has (T^) = {Ty) = 0. For large values of //. a product of 
neutron and proton BCS states is expected. Since this state is invariant under 
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the transformation exp(—i7r{A^,/2+Tz)), I enforce this symmetry, which entails 
A2 = 0. I furthermore assume that both gaps An = (—A^+iAy) / and Ap = 
( Aj. + iAy) I \/2 arc positive, as may always be achieved by a transformation 
of the form exp(— i(^Av + xTz))- When A is placed midway between the lowest 
and the highest efc, one then gets {A^ — 2Q and An = Ap = A for any value of 
due to the equidistant singlc-nucleon spectrum. To speed up the calculation, 
I keep A rather than G fixed with the variation of Then G varies in the 
case considered by less than .6 %. 

With quasinucleon annihilation operators = Y^kaMiMr^or + ^iMr^Xar) 
defined by 

[a^ , = , > 0, {ftj , a]} = b^^ , 
R^^Ho + A(P„ + + Pp + Ppt) + k{Tz)Tz - XA^ - fj,Tz, 

the RPA Routhian R is obtained by truncating to second order the boson 
expansion of R that results from making in the expressions for Hq, P, T and 
Ay the substitutions 

aja^ = bij + ---, ajaj = ^ blfiji, 

I 

where the boson annihilation operators b^j — —bj^ satisfy [bij,bi^] — and 
[bij,bi^] — 5ii5jm — Si^Sji- The normal mode annihilation operators P^ = 
T.i<j{(t>u,ijbij + '^u,iAj) are then given by 

[P,,P]=a;,P,, a;,>0, [P„ pt] = 5^^, 
and the ground state energy is [11] E = Eq + E2 with 




It may be noticed that when the terms with four quasinucleon annihilation 
operators or four quasinucleon creation operators are removed from the ex- 
pression for the Hamiltonian H so that the Bogulyubov vacuum becomes an 
cigcnstate of H (when the operators A^ and are replaced by their eigen- 
values in the relation H = R + XAy + A*Pz), then only the first term in the 
expression for E2 survives, and Eq + E2 becomes an exact expression for the 



3 




4, 



Fig. I. Eq — £Jo,T=0) E2 — ^2,r=o and /x/2 as functions of T. Eq — Eq^t=q is scaled 
by the factor 1/3. 

ground state energy. Therefore this expression derived from the boson expan- 
sion is expected to provide a reUable approximation even though most of the 
normal modes are not much different from two-quasinucleon excitations. 

The implications of the symmetries of this model for the normal modes are 
discussed by Ginoccio and Wcscner [12]. Two Goldstone modes result from 
the commutation relations [Av,-R] = [T^,/?] = 0. Furthermore, since [T+,i?] = 
/iT_|_ one normal mode has the frequency yU. Its annihilation operator is the 



linear boson part of T+/ y2{Tz), and it becomes a Goldstone mode in the limit 
// — > 0. The degree of freedom of this mode is the direction of the isospin. In 
particular the isospin quantum numbers Mt — T — (T^) may be assigned to 
the ground state of the Routhian R. The Goldstone modes contribute with 
the frequency zero to the expression for the second order energy E2. 

The parameters of the calculation are chosen so as to simulate the ^4 = 48 
isobaric chain: Q = 24, 77 = 2.1 McV, A = 1.7 MeV, k = 1.2 MeV. The 
result is shown in Figure 1. Eq — £'o,t=o depends essentially quadratically 
on T. It is in fact given in a very good approximation by the expression 
Eq — i?o,T=o = 1(^7 + k)T^ obtained for A = 0. The almost exactly quadratic 
T-dependence of Eq — -E'o,t=o is seen also indirectly from the hnearity oi fi — 
dEo/dT + 2(A/GfdG/dT, where the second term is negligible. E2 - E2,t=o 
shows a different behaviour. It rises for T ^ linearly with T and is in fact in 
this limit equal to /i/2. The linearity thus stems from the single term in the 
expression for E2 which represents the zero-point energy of the normal mode 
with the frequency /x, or, in other words, from the quantal fluctuation of the 
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isospin. 



The rest of the second order energy E2 is faily independent of T. This suggests 
that the other normal modes are to a large extent independent of the iso- 
rotational degree of freedom, or, stated otherwise, that the iso-rotation is 
highly collective. The deformation underlying this collectivity is in the pair 
field [13]. Thus with = the isovector A is perpendicular to the iso- 
rotational axis. So it breaks the iso-rotational invariance with repect to this 
axis, and a collective iso-rotation can arise. On the other hand the contribution 
to E2 from the non-collective modes varies from T = to T = 4 by almost 2 
MeV, so the RPA correlation energy should be taken into account in a detailed 
comparison of the results of Hartree-Fock-Bogolyubov calculations with the 
empirical masses. 

With Eo-Eo,T=o = l{v + i^)T^ and E2-E2,t=o ^ = IdE^/dT = \{r] + K)T 
wc have altogether E — Et=o — ^{t) + k)T{T + 1), that is, wc get the T- 
dependence of the symmetry energy found in the data. The Hartree-Fock- 
Bogolyubov energy expectation value includes only the first sum in the ex- 
pression for E2 [11]. This sum is cancelled to a large extent by parts of the 
second term in the second sum. For the contributions from the symmetry po- 
tential energy this cancellation is in fact exact. For A > the first sum 
in the expression for E2 is as a function of T analytic and even at T = 0. For 
A = it equals, however, ^T. Thus it produces in the absense of pairing a 
Wigner term, albeit with only n/ {rj + n) ~ 35 % of the full value. Although the 
forces there are different, this may explain the experience with Hartree-Fock- 
Bogolyubov and Hartree-Fock calculations mentioned in the introduction. 

A Wigner term corresponding to the term |T in the present model is actually 
the only one that may be derived from arguments like those in References 2- 
4, which are based on the form of the residual two-nucleon interaction. It is 
remarkable that with a deformation one gets also a term |T corresponding to 
the 'kinetic' part ^T^ of the symmetry energy. 

The symmetry potential energy of the present model differs from an 

interaction potential energy KY,i<jti ■ tj, where the index i or j labels the 
individual nucleons, by the term |«Av, which depends only on the number 
of valence nucleons. This interaction favours isoscalar nuclcon pairs. It 
is included in the Hamiltonian in order to get a realistic symmetry energy 
coefficient. It is in fact well known and discussed in detail by for example 
Bohr and Mottelson [4] that the kinetic term accounts for only a part of 
the empirical symmetry energy coefficient. Bohr and Mottelson consider a 
contribution to the single- nucleon potential energy corresponding to a two- 
nucleon potential of the form uti ■ t2. Thus the form of the symmetry potential 
energy of the present model is adopted from their discussion. This symmetry 
force is inessential, however, for the formation of the Wigner energy, which 
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results, as it was seen, from the collectivity of the iso-rotation due to the 
isobaric non-invariance of the isovector pair field. 



I am indebted to Stefan Prauendorf for drawing my attention to the issue of 
the Wigner energy and for discussions of the matter. 
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